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ocean model ere used to ccmpute isotherm displacemenzs at 
miant selected gtid peints in the North Pacific Ocean. The 
ocean model input parameters are climatological radiation, 
clouds, surface air tempémure, surface humidity and six- 
mourly synoptic winds taken from Pleet Numerical 
Oceancgrarphy Center analyses. BiCmErl ns tae baroclinic 
modes are comruted from the mean temperature profiles a* the 
@Cight selected mcdel grid points and these modes are then 
Ext in the time domain te the isctherm displacements. 

The fit shcews evidence of baroclinic free Rossby waves 
in the time demain at the selected grid points south of 4ON 
mmoeiundicates the opredcominance of the first barociinic medée. 
At the selected grid points north of 4ON, the results are 
Mreconclusive. This is due to the longer periods required at 
Bae Dorthern latitudes, tke neglect of salinity, and masking 
by a str¢eng seasonal signal caused by surface neating and 


cooling at higher latitudes. 
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I. INT ROD UCTION 


Sf 2am See we AE 81 ES SS ee 


QCceanic ceéeneral circulation models have bsen under 
d2velopment for the past decade anda half. EvVaiuation of 
these newly developed models was initially hindered by a 
lack cf adequate observational data from ocean crégicns. 
However, cver the past eight years high quality oc¢sanic data 
and data analysis studies from portions of the ocean have 
been made availakle to the large-scale ocean modeller. Of 


particular imecertance and interest in this stud has been 


tag 


the analyses cf temperature data resulting from the "North 
Pacific Experimrent® (NORPAX) begun during the last decade. 
Fased upon an analysis of TRANSPAC data (a@ sub-orogram 
SumechPaAXk Conrined to a région 30-50N and 170E - 150W), Kang 
and Magaard (1$80) deétergined that annual temperature fluc- 
Biats.ons in the area 30-40N, 160E - 130W consisted of firs= 
mode Rosscvy waves with random phasé traveling northwest wit! 
a wavelength cf 300 kn. HOG sone weicte te DdSs) 4Ot fhe Norch 
Pacific Océan, Price and Maqaard (1980) determined wnat 
first mode baroclinic Rossby waves with dominant period of 
S¢ven to t¢€n years existed in the region 40-50N. It was 
also determined that Rossty waves of annual pericd did not 


Mma crole in the latitude belt 4O-50N in the Nerth Pacific 
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Ocean. Tagadrome( i255), ih a2 Ppapsr discussing baroc 
Rossky wave energetics, summarizes these earlier studies by 
stating that first mode baroclinic Rossby waves dominate the 
internal temperature fluctuations in the region 20-30N in 
Mmemnorth Facific Ocean fcr a whole range of periods between 
five months and ten years. Between 30-40N only annual or 
néear-annual first mode bareclinic Rossby waves are evident. 
It is the resuits of these studies that have metivated this 
investigation, which is essentually an avaluation of whether 
a general circulation ocean model, driven by observed synop- 
eee winds, can simulate the above characteristics o£ the 
North Pacific Ocean Rossby wave field. Following a precs- 
dure outlined Ey Emery and Magaard (1976), this study exan- 
ines «he pessible existence of barociinic Rossby waves 
meer =sd frcm =amperature fields produced by the ocean 


model. 


aA 





II. MODEL DESCRIPTION AND ANALYSIS PROCEDUR 


A. MCDEL CESCRIFTION 

The mcdel used in this study isa threée-dimensiocnai 
primitive é€quaticn model developed at the Naval Postgraduates 
Schocl and described by Haney (1980). Atlas botzon and 
rigid iid are prescribed. The model is driven by seasonallv 
meeenmg (mencthiy climatolcgy) solar radiation, ciouds, sur- 
face air temperature, and surface relativ2 humidity. Tas 
surface winds utilized are taken from six-hourly synoptic 
analyses, 1969-1978, rade at the Flest Numerical 
Oceancgraphy Cénter, Monterey, California. Salinity is neg- 
lected. The mcdel domain shown in Figura 1 is a rectangular 
Basin extending from the Equator +o 65N and from 145E - 
125W. The herizontal grid spacing is approximately two 
degrees in latitude and three degrees in longitude. The 
model is 4000 m deep with 20 irregularly spaced ievels, 14 
of which ar¢ abeve 1000 x. The botndary conditions pre- 
Scribted at the bottom and sides are those of an enclosed 


thermally insulated basin. 


Vz 
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7 : , 
145°E 160°E 175, £ 170° WN 155 NW 140° W 125° 


Horizontal domain of temperature field used for analysis + 
shown by ¢nclosed regicn. Triangles mark geakd. ~ POAURt 
selected fer Rcssby wave anaiysis. 


Begare 1. Horizontal Domain of Ocean Model 


The mcdel also includes a dynamic adjustment mechanisn 
as descrited ty Adamec, et. al. (1981). Essencialily, he 
mechanism assures that if the local gradient Richardson nun- 
ber, Ri, falls below a critical value (0.25) then the temp- 
erature and tte velccity components at the two adjacen 


levels are mixed verticaliy so that heat and Momentum are 
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consé¢rved, +t2 mixing ratios for heat and momentum ere 
equal, and after mxing Ri equals the critical value. 


The model was run for a 30-year period, with «he 


- 


ih 
ct 


23s 


two decades considered asa "“spin-up"™ period and the thir 


fu 


decade, 1969-78, considered as a Simulation of actuéel condi- 


tions. The prirpcse of this thesis is to analyze the tampe 


'{ 


ature variability simulated by the model during this la 


in 
et 


decade. The cecal is to see if the variability is consistent 


With the dispersion relaticn fcr baroclinic Rossby waves. 


Be ANALYSIS FROCEDURE 
1. Geverning Equations fic 
The governing equations for baroclinic Rossby waves 
are derived from the frictionless momentum equations, the 


Igeee@ecstatic ctalance. and the ccntinuity equation for G 


pu) 


Pemeomeressible fluid in the follewing form; 


u, + C(V-V)u - fv + = Dy = Or 5 (la) 
a 
° + = =) -() (1b) 
v, + Sov jv = fu 5 Py, ) 
cee) | = Cee Cie) 
P, Eee = Oe (1d) 
PGT) o.= Oe. ire) 
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mmmere V iS the velocity vector, u, v and w are the velocity 
compcnents in the eastward (x), norvhward (y), and upward 


fpeeecirection, xrespectively, pis the density, p is the 


C7 


pressure, gis the acceleration due +o gravity, f is the 
: . : ‘ _ 3 re re 

Coriclis parameter, + is time and V= ( 3x" ay raz) ° The sub- 

Bememecs denote partial differentiation with respect to the 


subscripted variakle. 

iie 8 VGCrcticaty Squation is formed using egquaticns 
(1a) and (1b). Subsequent =) jojo sie eys else Oe Sane quasi- 
geostrophic approximation and linearization leads to the 


memetcwing Crm Of the vorticity equation: 


a Po pee : 
xX Buy BY Sees ee 05 | (2) 


where p' and wt are the perturbation pressure and vertical 
Merocity ccmecnent, B = df/fdy, p, is a constant reference 
density, and the other variables are as defined beéfors. 
Additional assumptions used in the above linearization pro- 
cess are that the lateral shear of tne pasic flowy is neg- 
lected and the basic state horizental velocity is a+ rest 
fewe., 2 = V = QO). From this assumption, equation (1c) and 


the Ecundary ccndition that w = 0 at z = 0 (sea Surface) it 


is 





eemescws that w= 0. Equation (2) contains two unknown vari- 
woees, p' and w*, therefore, an additional equation is 
required. 

Substituting from (1d) into (le) after multiplying 
by -g/p, and subsequent linearization yields the second 


equation, 


Z 
t t «= 
Peet ON wo = 0 Cea) 
Or 
a 
oS. 1 
w' = — Pot >? Csi) 
Po 
where N= - (J/p,) Pz is the Brunt-Vaisala frequency squared, 
assumed tc be a function cf zonly. Bliminating w' in (2) 


Msing (3b) yields, 
(pi, t+ pt), + ppt + (Sop = 0 C4) 
BOX Vii eo x NZ Potiz i 


Using the separation of variables «echnique the perturbation 


pressure is written, 


= 


p'(x,y,z,t) = Pi(z) py(x.y.t) , (5) 
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Where nm is the mede number. Substituting (5) into equaticns 


(4) and (3b) gives, 





27a a 
+ 
+ dit ode = oer t Pn’ x 4 ¢ (6) 
~ PP daz ye dz £203 +” i 2 
sb Po 2e 
and 
t 
medi = | Lo (7) 
Ne 2 (Bo. 
Po 1 
where Y*= ( a 8° d A is eo ti Stance wee 1 
= (ss2¢ gy2) ean mes BomeooPatatzioOn constant Wit! 
dimensions cf inverse velccity. Differentiating equation 
(6) with respect to z yislds, 
ac de Cae 
———- {— —"] + A eae U0 (8) 
2 GZ ce z 
dz N 


wgesc2cucing f£rem equation (7) into equation (8) yields the 


governing é€quaticn or vertical 2quation, 
-f é : 
down + a N?(2)w! = 770 (9a) 


eaoyec. tc boundary conditions, 


oe = «COO at z = O (Sea Surface) 
(9b) 


oa = 0 at z = -H (Sea Bottom) 


a 





meee ton (6) written in tke forn, 


V7 (Pade — AZ f(a de + BlPn, = 0 (10a) 


represents the Rcssby wave equation for which solutions of 


the fern 
ee = EXP[i(kx t ty - ot)} , (10b) 


are scugh+ where k and 2% are wavenumbers in the x and y 
direction, respectively and o is the frequency. Subs<es2ut— 
amg f£rcom (705) into (10a) yields a form of the dispersion 
memeetionshipr fcr baroclinic Ressby waves, viz., 
2 Z 
2 eee — OK + 2) 
d Se (11) 


n of- 


Equaticn sez (9) forms a Sturm-Liouville eigenvalue 
Peoplem for the vertical motion, where A, are the eigen- 
Meeeucs and w' the corresponding eigenfunctions. A solution 
is scught for a finite numkér of modes, m. These modes cor- 
respend tc the first m baroclinic Rossby wave modes. 

If N(z) were independent of depth at a fixed loca- 


tion and time, Pucmeccluctom@ foomew’) could be readily 
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determined analytically as SE LgOvOuetcl c PimMeetOns OF 
(n7z)/H (see for example Fedlosky (1979, pp. 356-360)), pro- 


vided 


J 
=} 


cuz) 


» 
2 
tt 
fi 
we 
ra 
i 
her 
ho 
WwW 


In certain instances an anélytical solution to equation ses 
(9) can be fcund provided N7z) 1S, cr approximates, a well 
defined furcticn (see for example Willmott and Mysak(1980) 
who solve the vertical prcklem with N4{z) as an exponential 
function Gite) Realistic distributions of N’are depen- 
dent upon depth and gqenerelly not well defined functions of 
@eeeetnerefcre the soluticn of equation set (9) is cbtained 


numerically. 


2. Determination of N* and The Isotherm Displ2cemen* 
Field 
The Brunt-Vaisala frequency squared is, as defined 
earlier, equal to -G/P,) Pp, ‘ Since the model neglects 


Salini«y, the simple stats relation, 


oli - a(T - i) (13) 


is used to convert the Brunt-Vaisala frequency squared frem 


eeeeruncticrn cf the vertical mean density gradient inte a 
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| 


mametsOrn Of the vertical mean temperature gradient, where 
fT, 1S a@ ccnstant reference temperature and a is the coeffi- 
cient of thermal expansion. Differentiating equation (13) 


Maeecespect tc Zand substituting into the definition of 


N? gives, 


n* = ag ips : (14) 


Equaticn (14) is used te calculate N? from the mean 
monthiy temperature fields output by the ocean model for the 
eet cen Years, 1969-78, cf the model run. Only che temper- 
ature field in the region 18.3 - 56.9N and 159.1F - 125W 
Shown as the enclosed box in Figure 1is used in the calcu- 
laticn of nt After completion of the numerical differentia- 
mon (ri~ght side of egquaticn (14)), a ten-year time mean, N2 
at each grid point and at 19 vertical levels (mid-levels of 
the Z0 model levels) ‘is calculated. Bana y Nis averaged 
Over 28 seast-west grid points, giving a N? versus depen. pre= 
file for each latitude (2¢0 north-seuth grid points). Justi- 


fication fcr the east-west averaging is based upon 1i++132 


wD 


variability of N?with longitude over the region at a partic- 


ular vertical ievel. The variability of N2 with latitude, on 
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the cther hand, womcllcnmienat gOrch-south “averaging ic not 
justified. Price and Magaard (1980) also found that mean N? 
profiles varied much more with latitude than with lengitude, 
and thus the numerical mcdel daéta are consistent with the 
corresponding data taken from the Nerth Pacific Ocsan. Pro- 
files of (N2 y/ , averaged over longitude, ar2 shown in 
Figures 2 to 5 for the latitudes of the grid points selected 
in this investigation (triangular symbols on Figure 1). The 
solid curves cn each figure are drawn to the values of ("2/2 
at 19 vertical levels determined Ey agquation (14). The 
dashed curves are drawn to values of (N2)/2 Galewlat2daa= 
every five meters from z = 0 to z = ~4000 m uSing interpola- 
mon . The interpolation is done numerically by means of 
Aitken's scheme of Lagrange interpolation described by 
Hildebrand (1956, pp. 49-50) using the values of (v2) at the 
19 vertical levels and an 4a@ssigned value of (n2 7 at -4000 on 
Gemeeakular input. At latitude 48.8N and 54.8N the assigned 
bot~tem value for (N2)/? Was given the identical value calcu- 
lated at level 19 (2830 m). At the other salectec latitudes 
*he assigned Ecttom value was determined by linear extrapo- 


laticn using the levels below 1000 mn (levels 15 to 19). 
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[MiwNEErOfrLiles Sh¢wa pronounced variability with 
latitude. Belicw approximately 900 a, (n2 y/? is small (less 
mmanm 0.001 1/sec) and nearly constant (model simulated 
isothermal ccndition in the deep ocean). At latitudes north 
or 47N this constancy is evident at depths béiow 150 m. The 


maximum value cf (N2)/? 


(Taximum vertical temperature gradi- 
mee pocccurs at level 2 (16.5 m) at latitudes south of 22N 
and Fetween latitude 4ON and SON. At the other la*itudes 
the maximum cccurs near the surface (5 m). A secondary max- 
Sum occurs between latitudes 30N and 45N at a depth between 
200 m and 400 @ (inferred lccaticn of the mean thermocline). 

The secord field te be ccmputed using the model sin- 
ulated mean mcnthly temperature fields is the isotherm dis- 
placement. The isotherm displacement, D(Zet), iS tne 
displacement of @ monthly mean isotherm from its long-tern 
average level. It is defined mathematically as -T'/ (T), 
where T' is the temperature fluctuation, T' = T - T. Vand 
(T), are the long-term (ten-year) averaged temperature and 
vertical temperature gradiert, respectively. In the numeri- 
Cal calculaticn of D(z,t) the assumption is made that tne 


tempereture field is linear in z between mode2l levels. The 


isotherm displacement versus time curves at 19 vertical 
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levels Efrem tke ten-year mcdel Simulation are shown in 
Mmeagurss 6, 7, and 8 for the grid pcint at 30.5N 139.1%. 
Ppmveectic Came ts —deranitlon, Di{z,t) positive, ir 
general, imciies a temperature cclder than the long-term 
mean. D(z,t) negative implies a temperature warmer than the 
long-term mean. The isctherm displacement curves, as 
depicted in Figures 6 to 8, reveal a pronounced seasonal 
Signal (annual period) at levels above 71.5 nm. At 100 o 
(level 6) and Frelow the seascnal signal becomes less well 
defined and lcsses definition with depth. The seasonal sig- 
nal at the upper levels is primarily the model sinulated 
mean monthly temperature fields response to the seasonally 
varying solar radiation prcevided as model input. Below 100 
m (Figures 7 and 8) a low frequency (long period) signal is 
evident. The isotherm displacement curves at the other 
Selected grid feints (not shown), in general, show similar 


Characteristics cf an upper level amnual cycle and low fre- 


quency signal at lower levels. The Significant features, 
Meeocecding frcem south to north at longitude 155.94, are 
eivat , 1) the depth at which the seasonal signal remains 
distinguishable increas2s reaching 1760 m (ievel 17) one 


54.8N, 2) che amplitude cf the isotherm displacemen* at 


Ze 





mid-l¢evels increases from a maximum of approximately 50 m at 
meomon tO 5000 m at 5S4.8N. This increase in amplitude is a 
result of the néar iscthermai conditions evident below 150 o 
Meethe nerthern latitudes (i.e. OB) is small thereby 


increasing the value cf D(z,t)). 


(a) 


Equaticn set (9) is solved numerically using 23n 
Interactive cCrdinary Differential Equation (IODE) solver 
which uses a fcurth-order Runge-Kutta integration technique. 
The eigenvalues, eee So -eeecCunGmoye Mans Of —zhe "c=reor 
methed" descriked by Preisendorfer and Gonzalez(1973). Ai) 


Meieeying the "c-rcot methcd" equaticn (9a) is separated int 


mm@omtaest—-Ccrder differential equations, 


I 
= 


(15a) 


Oo. 
N 


: 


x2 N2czdwt (15b) 


n 


Se 
N|o 
t 
i 


A triai value is choosen for A, and aquations (15a) and 
Mpeo) are integrated using IODE starting at ~he sea botton, 
Z = -H (4000 m), with boundary vaives of w' = 0 and 9= 1. 


n 


As the integration proceeds the zeros of the function 
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function C(Zy = w'(zZ) - 9(2)2(z*) are noted. econe 
above expressicn 2(Z') equals w'(z"')/g(z') where z' is the 
value of Zz at the bottom kcundary which changes as =~he inte- 
Geeeen continues. Since the boundary condition states thac 
Weep = w' (zt) = 0, then Z(z) must equal zero and the c 
function reduces to c(z) = w'(z). The eigenvalue problem 
thus reduces tc finding values of z where w'= 0. 

meE example, for the (N22 DEOtLv a Maews0s5N and a 
Meee Value CL A= 2.5 sec/m, c(z) = wt (z) = 9)when z= 
Seem, —~1-:C&m, -607m, -326m, and ~112n. Niaeey PSG Xo DY Oye. | 


maeym, noting when w'= 0, and tracing out the points result 


7) 


in five curves (eigencurves) as shown in Figure 9. Jee 


iD 


first five eigenvalues which satisfy the boundary condi- 


tions, equeticns (9b), are found by the intersection of the 


cf 


eigencurves with the z = 0 line. The firs: five eigenvalues 
Mechs E€xample are shown to the right of the ordinate in 
Figure 9; namely, A;,= 0.4735 sec/r, A, = 0.8934 sec/n, 

eee 1.3380 sec/a, Ag = 1.7814 sec/n, Ag = 2.1823 sec/n. 


Meer t2caticn for limiting ths calculation cf eigenvalues to 


the first five mcdes will be made apparent in the next chap- 


Mm 


Bec. An advantage in finding the ¢igenvalues in this manner 


xed. bb 


ity 


is that the depth to the sea bottom need not be 


a2 





mmemaepth to the brottcm were, for example, 3 kw rather than 
4 kn, the cigenvalues satisfying the same given boundary 
@emaitzons could be oktained by noting the valué of at the 
intersecticn cf the Eigencurves and the z= -1000 am line. 
The recirrocai ¢igenvalues, Shea lLen) COREeSDONdG €o the 
internal phase speeds in the absence of rotation associated 
With each kEarcclinic mode, are given in Table I for tha 
latitudes of each selected grid pcint. The tabular values 
Of c, show a decrease with latituce and with increasing mcde 
number. 

Using expendartle bathytkermograph (XBT) data ctaken at 
weather staticn November (30N, 140W) during July 1966 to 
December 1569, Emery and Magaard (1976) obtained c, values 
Of Cc, = 2.99 m/sec, Co = 1.41 m/sec, cy, = 9.92 m/sec, Cy = 
0.72 m/sec, and co = 0.51 m/sec. Thes2 values ar2 slightly 
higher than the values given in Tabl2 I for latitude 30.5N 
Measatude cf amcdel grid feint closest to the liatitude cf 
weather staticn November). The difference can be attributed 
to the deeper depth at weather station November (4755 mn) and 
#0 larger values for (eae the vaines for N*calculated 
from «the modél simulated temperature field is smaller than 


observed). 
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Having determined the first five eigenvalues by the 


N"o-rcot method", «ne next step is t9 calcui 


respcendina 


eigenfunctions, 
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an 


; ise 0a 
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TABLE. 4 


Reciprccal Eigenvalues (m/sec) for 


Acdes 1 to 


a 


easeETODE MODs NUMBER 
—— 1 ol ae Ee 4 > 
18 .3N DIST 1.2442 One 50'5 0.6274 9.5049 
24 .4N Cece | We Uso 0.8202 076079 0.4922 
Sis ON 2.1119 Nev los 0.7474 0.5614 9.4582 
36.6N N79 37 0.9774 0.6245 0.4719 Ole Sieltels 
42.7N 6 eVAlS, €575 10 0.4736 8 Silas 7a, 0235077 
48.8N ol BSNS) 0.4831 0.3294 Omezo2 9 0.2264 
54. &8N 0255 11 UecoS Ue los 0.1722 DasS © 
@elgenfuncticns are then nermalized, following Emery and 


Magaard (1976), as 


wi ' Cz) 


feats? = (16) 





Oo ? 
( / WCAG Ea) erat 


-H 
where n represents the mcde number. 
Tepresents a Sturm-Licuviile system, this 
Seencgonal with respect to the wezght function N 
eee ve>fical cclumn z = -H to z = 0, vViz., 


Oo 


| N*(z)w'Cz)w,"(z)dz =O. n #Mm alge, 


-H 


Shs 





Therefore, ncrmalization of the eigenfunctions results in 


the crthonornal S62 05 14) ae ee 


oe. 
a edz )y.(2)y_(2) dz = 6, (18) 
-H 

where 6, is the Kronecker delta function. A more détailed 


MMecussicen Concerning the orthogonaiity condition for ths 


elgenfunctions ls given by Willmott and Mysak (1980). The 


normalized sigenfunctions for the selectad latitudes are 
shown in Figures 10 through 16. The characteristics of the 


Smeemcunctions ares related to the profile of nN? (Zia. At 


ae) 


latitudes south of 4ON (Figures 10 to 13) th 


iD 


BOda evo. nTs 
of the normalized eigenfunctions are concentrated near the 
upper half of the column in response td the model simulated 
upper ocean N? maximun (Maximum temperature gradient). At 
the northern latitudes (Figures 15 and 16), the normalized 
6igerfunctions are very nearly sinusoidal, particularly at 
the higher omcdes. This corresponds to the near constant 
Beotils of (2? (Figures 4 and 5) below 150 ao. 

The analysis procedure is concluded by fitting the 
normalized eigenfunctions to the simulated isotherm dis- 


placement, 0D(Z,%~)- nic 296. dOmein the <ime domain dy 
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making use ct the cempiezteness of che eigenfuncticns <<o 
equate D(Z,t) to the sum of the normalized eigenfuncticns 


Seeeipiited by an amplitude coefficient A, viz., 


Diez. t)- = y BCE) a) a €799 
n=1 


Mmemecertficients, A,(t), are readily obtained by making use 
Seeene Orthcacrmality condition. Multiplying each side cf 
Peation (19) bv N*Y,,end then integrating from z = -H to z= 


QO resuits in 


oF. 


e) 
/ N? (zy (z)D(z ,t dz : sD [A_(t) | N*(2)¥_(2)¥_( 29d] (20) 
H 


= n=l 4 


Mops cat21on cf tke othogonality condition (18) to equation 


(20) yields 


‘o) 
: 2 : ; 
A, t) = | N (zy (2) D(z ,t)dz : Ca), 
H 


Bejgeeden (21) is vsed to cempute A(t) for the tarst 
meve modes for each sslected grid point. The daviation from 


mime mcdal fit, as defined Ey Emery and Magaard (1976) ,eqn(22), 





SemassoO Caiculated. The deviaticn term quantitatively snows 
how well the model simulated isotherm displacemen= fluctua- 
tions are fit by the first five baroclinic modes. A value 
of DEV(t) = 0 implies that the isotherm displacement fluctu- 


mewens can be fully explained by the first five baroclinic 


modes. 
ed 5 
: N2(z)([DCz,t) - > A (ty (2) )? dz 
n= 1 
re ee ) 
@) 
/ N Cz) Deeds 
-H 
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ieee RES Ui: 


lun 


Wemcocriscients, A,(t), and the deviation, DEV(+t), cal- 
culated at each selected grid point are shown in figures 17 
po 24. At lower latitudes (Figures 17 to 19 and Figure 24) 
the amplitude cf coefficient A, is generally larger than the 
S@eeercients Ao tc A. “ndicating the importance of the first 
mode at «che scuthern latitudes. At latitude 24.4N (Figure 
eyeeae Seascnal signal is evident for coefficients A, and 4.. 
The seasonal signal of the higher modes becomes more appar- 
ent with incr¢easing latitude indicating that the seasonally 
varying mecdel input parameters appear to be reflected in the 
higher modes at higher latitudes. At the northern most gri 
Beemc, O4.EN (Figure 23) the overall amplitude increases an 
order of magnitude. Additionally, the pattern of each ccef- 
Mmement afrears neariy identical with only siight differ- 
ences in magnitude. 

The deviation term at latitudes 18.3N and 24.4N remains 
generally belew 50% =during most of the ten-year model run, 
althcugh scne time periods are higher. Nerth of 24.4N there 


is a general increase cver the ten-year period in the 
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Meegare 21. A, (t) and DEV(t) at 42.7N 155.98 
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Figure 22. aA,(+) and DEV(t) at 48.8N 155.94 
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Figure 23. A,(t) and DEV(t) at 54.8N 155.9W 
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Figure 24. A,(t) and DEV(*t) at 30.5N 139.18 
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@@viation curve with the excepticn of the two northern acst 


th 


Geaa Points (Figures 22 and 23). AS a comparison, the first 


Lj 


five coefficients and the deviation term computed by 2mery 


and Magaard (1976) using 48 months of data fer weather sta- 


tion Novemkeéer are shown in Figure 25. The results shown in 
imeemcres 24 (model grid rcroint closest £90 weather station 
November) and 25 compare reasonably well. Both show the 
meadest amplitude of coefficeint A,. The 48-mceath mean 


deviation ccmputed for weather station November was 32% 
while the ten-year mean deviation is 38.2% for model grid 
Been. 30.5N 139. 1f8. 

To examine the relative contribution made by each dbaro- 
clinic mode tc fittizrg tke simulated isotherm displacement 
Mmictuacticns (eguation(t9)}), the ten-year mean deviaticn 
(2@quation (22)) was calculated for each selected grid point 
adding one mode ata time. The results ac? given in Tacle 
it. Column 1 (under modes) gives the ten-year mean 
Gevaation when only the first baroclinic mod¢ 1s iit to the 
Msotherm displacement fluctuations. Column 2 gives’ tie 
ten-year méan deviaticn when modes 1 and 2 are fitted to the 
isotherm displacement fluctuations and so on. The tabular 


values show ‘the importance of mode 1 at the selected grid 
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Figure 25. an (t), and DEV (t) for weather station November 
(from Emery and Magaard, 1976) 

memes) sOuth cf 30.5N (1.¢€., over 40% of the deviation is 

eumeetunted for ty the first baroclinic mods alons). The cca- 

Mme On made by the fifth barcclinic mode in reducing the 

mean deviation is less than 6% at the thr2e southern mcsz 


Mec points, and, except fcr the northern most grid point, 


’ 
) 


, ’ 
cl 
U 


(1) 


is less than 1C% elsewheré. Based upon these results, 
apparent that adding higher baroclinic modes (modes gr2ater 


meen 5) weculd adG@ little te reducing the deviation term conm- 


pared to mcde number cne. 
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TABLE If 


Mean Teviation (%) at Selected Grid Points 


GRID FOINT MODES 
1 1+ 2 i Z+s 1+ 24344 1+#2+3+4+445 
£8.3N 155.9W 2653 46.0 36.4 28.6 PAS as 
merean 155.98 Se 45.9 B20 B5e0 31.6 
omen 155.9N 57.4 $4.2 SUAS 46.9 41.0 
Soe ON 155.9W See Sheds scree. 45.3 Eisele 
BeeeN 155.958 Gr. 2 Eos 1 61027 54.2 Loe2 
Woe oN 155.9 56. | 34.5 2826 24.1 Bo 
Secon 155.98 oC ee str as 43.0 2626 Los 0 
Soon 139.100 oa (6 SOe5 30 As, 41.1 38.2 
The mean deviation (column 5 of Tabie fM[T) generally 


increases with latitude ur to 42.7N but then decreases at 
the two highest latitudes. Adds sonedny, the GONntribu< on 
Of medes 2 tc 5 in reducing «he mean deviation is signifi- 
cant at the ncerthern most foints. To 2xamine the relevance 
of these resuits, the cutoff period (minimum pericd) for the 
propagaticn ci free Rossby waves of mode n (see Fmery and 
Magaard (1976)) given by T = (47f4A,)/B was evaluated. The 
results are shcewn in Figure 26. The dotted line Gelinéates 


the annual period and the dashed line the five year pericd. 
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Figure 26. Cutoff Period versus Latitude 
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The curves cf Figure 26 shcw the following results: 


*® Pased upen the theoretical cu 
raroclinic free Rossby waves n 
Only mode 1 can cccur. 


f= period annual period 
€xnet NOren Cf 25N but 


e Paroclinic free Rossby waves at the annual period can- 
not exist north of approximately 36N (recall Magaard 
(ges) Seated that first mode baroclinic Rare! neues 

TO 4 Lo) ei 


were eae evident between latitudes 30 
pieces NCrt aClELC Ccean):.< 


py 
< 
Dp 
09) 
fw 
ct 
a 
Se 


e Barociinic free Ros 
@ainect Exist north c 


The results given above are significant. Since the frze 
ROssky wave analysis was conducted using only ten years of 
model simulated temperature data, free Rossby waves of peri- 
ods greater than approximately five years can not pe ade- 
quately resolved. Therefcre, the results obtained for model 
Meroe coin® S4.EN 155.9W lcsses significance. Likewise, at 
Mma point 43.6N 155.9W the theoretical cucofft period curves 
show that oniy firs: mode t-aroclinic Rossby waves at periods 
greater than approximately three years can ¢xist. Ts 


Bere, the low deviation at latitude 48.8N resulting prima- 


Mey trom modes 2 to 5 is not necessarily due to baroclinic 
Rossty waves. The opposite may be said of the seisccted grid 
Beents a= latitudes 13.3N *o 36. 6N. At these latitudes th2 


feest baroclinic mode is dceminant and the cutoff veriod is 
4.5 months at 18.3N and 12.8 months at latitude 36.6N. At 


Meareitude 42.7N the cutoff period is 21.4 months. HOwWEVEr, 





only the first and second modes can exist at periods less 
than five years. The mean deviation is high when only the 
First and seccnd modes are fitted to the isotherm displace- 
ment fluctuations (96.2% and 88.1% in columns 1 and 2, 
respectively, of Table II). Thus at model grid point 42.7N 
fetes coW che macdel Simulatsi isotherm displacement fluctua- 
tions cannot be explained ty free baroclinic Rossby waves in 
the time dcmain. 

Thus far the examinaticn of barociinic fr2e Rossby wave 
Signature in the model sinulated temperature field (or iso- 
therm displacement) has teen limited to the time domain 
oa ly. In order to conduct an analysis of baroclinic free 
Rosskty wave signature in the space domain requires carrying 
Bemea CLOSS~—-sprectral £3t cf the x and y direction wave xzun- 
bers at «he apprepriate frequencies. In the context of “his 
Study a cursory #xamination of the wave number vectors will 
be made. The iccus of wave number vectors for mode 1 at 
apprepriate frequencies fcr latitudes 18.3N and 36.6N are 
Shewn in Figures 27 and 28. The locus of wave numbers was 
determined using equation (11) for three different frequen- 


cies. 
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Figure 27. Lecus of Wave Number Vectors for Mode 1 at 18.3N 
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Figure 28. Locus of Wave Number Vectors for Acde 71 at 36.6N 





The spacing retween model grid points utilized in this 
SieayetS Cr tbe ecrder of 226 km in the north/south direction 
and 251 km (at 36N) in the easty/west direction. Tne shcr- 
test wavelength resolvable by the modal is equal to twice 
the spacing ketween grid foincs, which in this case is on 
the crder cf approximately 500 km. Therefore, Rossby waves 
resolvable by the model must have wavelengths greater than 
Or egual to the 2ad value, where Ad is the grid spacing. 
The wave numbers in the x ard y direction which correspond 
to a wavelenath which is greater than tne 24d waveiength are 
located withir the dotted rectangle in Figures 27 and 28. 
Stated ancther way, wave numbers cutside of the boxed-in 
regicn corresrcend to wavelengths too small to be resolved by 
the mcdel. 

The innexz circle in Figure 27 is the locus of wave nuna- 
bers found when a frequency corresponding *0 a six-monta 
period is uséd. Phew aetle cit cle and SuUTersezrcle are for 
frequencies ccrresponding to periods of nine months and 12 
months, respectively. Note that due to the different scale 
of eacn axis the circles appear as ellipses. Pat Teens. or 
all three circles fall within the enclosed regicn indicating 


that wavelengths greater than the 24d wavelength can occur 
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at the aforementicned freguenciss (periods). The circles in 
Figure 28 are fcr frequencies corresponding to periods of 18 
months (inner circle), 24 months (middle circle) and 30 
months (cuter circls), respectively. It is noted that at 
36.6N only waves having a period of about 30 menths or 
greater have wavelengths long encugh to be resolvable by thea 


model. It fcellcws that waves of 4éad or B8Ad would have even 


greater periods. Therefore, increasing latitudé requires an 
increase in zke Ressby wave period to allow for wavelengths 
long enough tc be resolved by the model grid spacing. Addi- 
tionally, using equation (11), it may be snown that for 
wavelengths that are resolved by the model grid, k7 + Q%& A2 £2 
and thus the period 27 /g must be larger taan the cutcft 
beticd. At high latitudes such long periods would not he 
well resolved ty the ten-year simulation. Thus, even though 
low deviaticns result at high latitudes (column 5 of Table 
EI) the isotherm displacement variability cannot be attrib- 


uted to actual Rossby waves with any degree of certainty. 
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ive CONCTUSICNS AND RECOMMENDATIONS 


Pased upcn the preceeding analysis «he following results 
are summarized: 


roclinic Rossby wave signature 
LsOtherm displacement fluctua- 


e On the average, fre b 
ula? 
! eae Puiemdcnesn FOG es tneuce |seted 


in the. mcdel sin 
tions is evident i 
Seed ECCintS south o 


}-- 


e At the selécted arid poz MicSmeSOUch On MYON slows eaLes 
Mrociinic node shows predominance. 


e Based upon the locus of wave number ellipses, Ro 
wAVEeS at the Ee pbee Tiate periods are long enough + 
resclvable w: he model grid spacing. 


® Fased upen the deviation term, certain times within the 
ten-year model run period indicate that the first five 
baroclinic modes do not adequately fit the isothera 
d-:splacement fluctuations (i.é€., the deviation is high, 
Bee rore the isotherm displacement Sees os can- 
not fully S“puained by baroclinic fre Rossby 
eee) 

e At selected grid pcints north of 40N the analysis 
results are inconclusive. 

e There is eee = bly gccd Sue reement between the analysis 
results at me ofealie! [exe ie DUee te loos hHee and cheats con — 
ducted by paery a aay Magaard (1976) for weather station 
November (3CN " 


The results stated in point number five above, are 2 product 
Smeceveral factors. The model run tim? (ten y¢ars) sets a 
constraint upen the periods that may be adequately resolved 
by the analysis method. The cutoff periods of tne selected 
Gad “Points ncrth cf 40N gradually begin to approach a 


five-year period, and wavelengths resoived by the model grid 
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Mesult im pericds considerably greater than che cutoff (nin- 


imum) period. Periods greater than fiv2 years can not be 
adequately resclved by the analysis method. A longer run 
time is reguired. Secondly, since the ocean model is a 


"fresh" water model the effect is to dimenish «he nagqnitude 
of the Brunt—Vdisdla frequency, Bat cules Vcc ioh= 
latitudes. This results in eigenvalues that are large conm- 


temeaeco "real" ccean values (Ais inversely proportional to 


N) . The larger 2igenvalues then, in turn, increase the 
value of the cutorf pericd more than what cccur when cnly 
the Coriolis parameter is considered. Poe ee Stem onal 


Signal, resulting from the seasonally varying model input 
parameters, dcminates the isotherm displacement fluctuations 
OVGI a greater portion cf the water column az the northern 
selected grid feints. Thus, the seasonal signal, which is 
primarily a thermodynamic effect, contaminates the first 
five modes. 

Eased upon these results several recommendations for 
further research may be set forth: 


e The seéscnal Signal resulting from thermodynamic 
effects snould be remeved. This could be achieved by 
€ither neglecting the model layers above the main ther- 
mccline cr ky applying a filter technique to remove the 
the ce seasonal signal. |The former suggeszion 
aoecmoomeuteeactscal at she high latitude model grid 

Poln=s, hOW2VET ¢ due to the penetration of the seasonat 

Signal dcwn to the bettom of t water column. 
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2 es, ScUCy Set tmepmcean curtbens equal to zero. Further 
studies sheculd include the mean current in ~*+he analy- 
Pes ee Dat tlcUulansly 2 che ©Fegion of the Caliltcrnaia cur- 
Meme, ana Necth Pacific Current. A study by Kang and 
Poadacrd (i279) has Shown that the North Pacific Current 
mised SigPitican= impact on the first baroclinic shear 
mode. 

e Since tke model neglects salinity furthe res2arci 


pares jpg 

ememeea study the Eéeasability cl incorporating salinity. 
This would frovide fcr more realistic density profiles 
and. improve the analysis results at =n 
latitudes. 


e The mech 
Eaclt ic 
ysis he 


anisms which produce Rossby waves in the North 
CGCeameiS Stita ne= Euily understood. Ths anal- 
Shown that certain tims Deere within £h¢ 
model ru Pena “clei lOcelacscude  4ON. cannot be 
explained fully bY Pie syste tive baroclinic modes. 
apeascently some other factor is involved. Further 
Research could be done to isolate any peculiarities ia 
tne the syncptic wind field used as model input Since 
Peewee ciX-hourly Syneptic Wind was the ocnly varyin 


rf} U) 
Ha | 


erve) 


tarametsr from one year to the next. The study shoul 
also te extended to investigate the source of Rossby 
waves (i.¢., whether they tesult fron large-amplitude 


Sweetie ns Cr the wind Stress curl near the Ssastern 
Ecundary as hypothesized by White and Saur (1981) of 
whether scme Other mechanism is involved). 


Sent. furtehsr studies shculd includ® comparison to apprto- 
pane <e cbservational studies when availadle, and 
hereky include a model evaluation. 


Thas study has keen a first step in analyzing the low fre- 
quency dynamics of the Haney (1980) ocean model below the 
mixed layer. A complete understanding of the dynamics of 
the large-scale variability of the ocean temperature field 
below the mixed layer (thermociilne) is the primary goa. 
This understanding will in turn lead to improvements in the 
Ocean model and its diagnecstic and prognostic capabilities, 


the ultimate gcal. 
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